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ABSTRACT 

Peculiar velocities of clusters of galaxies can be measured by studying the fluc- 
tuations in the cosmic microwave background (CMB) generated by the scattering 
of the microwave photons by the hot X-ray emitting gas inside clusters. While 
for individual clusters such measurements result in large errors, a large statistical 
sample of clusters allows one to study cumulative quantities dominated by the 
overall bulk flow of the sample with the statistical errors integrating down. We 
present results from such a measurement using the largest all-sky X-ray cluster 
catalog combined to date and the 3-year WMAP CMB data. We find a strong and 
coherent bulk flow on scales out to at least 1 300/;.~^Mpc, the limit of our catalog. 
This flow is difficult to explain by gravitational evolution within the framework 
of the concordance ACDM model and may be indicative of the tilt exerted across 
the entire current horizon by far-away pre-inflationary inhomogeneities. 

Subject headings: cosmology: observations - cosmic microwave background - early 
Universe - large-scale structure of universe 

In the gravitational-instability picture peculiar velocities probe directly the peculiar 
gravitational potential [e.g. Kashlinsky & Jones 1991, Strauss & Willick 1995]. Inflation- 
based theories, such as the concordance ACDM model, predict that, on scales outside 
the horizon during the radiation-dominated era, the peculiar density field remained in the 
Harrison-Zeldovich regime set during infiationary epoch and on these scales, the peculiar 
bulk velocity due to gravitational instability should decrease as K-ms oc and be quite 
small. Peculiar velocities can be obtained from the kinematic SZ (KSZ) effect on the CMB 
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photons by the hot gas in clusters of galaxies [e.g. Birkinshaw 1999]. For each cluster the 
KSZ term is small, but measuring a quantity derived from CMB data for a sizeable ensemble 
of many clusters moving at a coherent bulk flow can, however, overcome this limitation. As 
proposed by Kashlinsky & Atrio-Barandela (2000, KA-B), such a measurement will be dom- 
inated by the bulk-flow KSZ component with other contributions integrating down. This 
quantity, the dipole of the cumulative CMB temperature field evaluated at cluster positions, is 
used in this investigation of the 3-year WMAP data together with the largest X-ray selected 
sample of clusters to date to obtain the best measurement yet of bulk flows out to scales of 
^300/i~^Mpc. Technical details of the analysis are given in the companion paper (Kashlin- 
sky et al 2008 - KA-BKE). Our findings imply that the Universe has a surprisingly coherent 
bulk motion out to at least ^ 300/i~^Mpc and with a fairly high amplitude of 600-1000 
km/sec, necessary to produce the measured amplitude of the dipole signal of ~2-3/LiK. Such 
a motion is difficult to account for by gravitational instability within the framework of the 
standard concordance ACDM cosmology but could be explained by the gravitational pull of 
pre-infiationary remnants located well outside the present-day horizon. 



1. Method and data preparation 

If a cluster at angular position y has the line-of-sight velocity v with respect to the CMB, 
the CMB fluctuation caused by the SZ effect at frequency u at this position will be Si,{y) = 
STSz{y)G{u) + SKsz{y)H{u), with 5TSz='rTx/Te,ann and SKsz=rv/c. Here G{u) ~ -1.85 to 
— 1.25 and H{iy) ~ 1 over the WMAP frequencies, r is the projected optical depth due to 
Compton scattering, Tx is the temperature of the intra-cluster gas, and /cB^e,ann=511 keV. 
Averaged over many isotropically distributed clusters moving at a significant bulk velocity 
with respect to the CMB, the dipole from the kinematic term will dominate, allowing a 
measurement of Huik- Thus KA-B suggested measuring the dipole component of Si,{y). 

We use a normalized notation for the dipole power Ci, such that a coherent motion 
at velocity l^uik leads to Ci^kin = ^CMB(''')^Kuik/c^) where Tcmb = 2.725K. For reference, 
•\/ C'i,kin — l((''')/10~^)(T4)uik/100km/sec) /jK. When computed from the total of A^ci posi- 
tions, the dipole will also have positive contributions from 1) instrument noise, 2) the thermal 
SZ (TSZ) component, 3) the cosmological CMB fluctuation component arising from the last- 
scattering surface, and 4) the various foreground components within the WMAP frequency 
range. The last of these contributions can be signiflcant at the lowest WMAP frequencies 
(channels K & Ka) and, hence, we restrict this analysis to the WMAP Channels Q, V & 
W which have negligible foreground contributions. The contributions to the dipole from the 
above terms can be estimated as {S^{y) cos 9) at the A^ci different cluster locations with polar 
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angle 9. For A^ci > 1 the dipole of 6^, becomes aim - af^ + «?m^ + «?m ^ + ^7#^- Here af^^ 
is the residual dipole produced at the cluster locations by the primordial CMB anisotropies. 
The dipole power is Ci = Ylm=~i kimP- The notation for aim is such that m=0, 1, — 1 
correspond to the (x, y, z) components, with z running perpendicular to the Galactic plane 
towards the NGP, and (x, y) being the Galactic plane with the x-axis passing through the 
Galactic center. This dipole signal should not be confused with the "global CMB dipole" 
that arises from our local motion relative to the CMB. The kinematic signal investigated 
here does not contribute significantly to the "global CMB dipole" arising from only a small 
number of pixels. When the latter is subtracted from the original CMB maps, only a small 
fraction, ~ (Ai'd/A^totai) ~ 10~^, of the kinematic signal Ci,kin is removed. 

The TSZ dipole for a random cluster distribution is aj^ ~ ((rTx)/Te,ann)^c/^^ 
creasing with increasing Nc\. This decrease could be altered if clusters are not distributed 
randomly and there may be so me cross-talk between the monopole and dipole terms espe- 



cially for small/sparse samples (jWatkins fc Feldman 19951 ). but the value of the TSZ dipole 



will be estimated directly from the maps as discussed below and in greater detail in Kashlin- 
sky et al (2008 - KA-BKE). The residual CMB dipole, Ci^cmb, will exceed 0"cmb/-^c1 because 
the intrinsic cosmological CMB anisotropies are correlated. On the smallest angular scales 
in the WMAP data ctcmb — 80/iK, so these anisotropies could be seen as the largest dipole 
noise source. However, because the power spectrum of the underlying CMB anisotropies is 
accurately known, this component can be removed with a filter described next. 

To remove the cosmological CMB anisotropies we filtered each channel maps separately 
with the Wiener filter as follows. With the known power spectrum of the cosmological CMB 
fluctuations, C^*-^^^, a filter in £-space which minimizes {{5T — (^instrument noise) ^) in the 
presence of instrument noise is given by Fi = (Ci — C^^^^) / Ci, with Cg being the measured 
power spectrum of each map. Convolving the maps with Fi minimizes the contribution 
of the cosmological CMB to the dipole. The maps for each of the eight WMAP channels 
were thus processed a s follows: 1) for Cf *"^^ we adopted the best-fit cosmological model 



for the WMAP data (IHinshaw et al 20071 ) available from http:/ /lambda.gsfc. nasa .gov; 2) 
each map was decomposed into multipoles, aim, using HEALPix (iGorski et al 20051 ): 3) the 
power spectrum of each map, Ci, was then computed and Fi constructed; 4) the aim maps 
were multiplied by Fi and Fourier-transformed back into the angular space {0,(f)). We then 
removed the intrinsic dipole, quadrupole and octupole. The filtering affects the effective 
value of T for each cluster and we calculate this amount later. 

Here we use an all-sky cluster sample created by c ombining the ROSAT-ESO Flux 



Limited X-ray catalog (REFLEX) (jBohringer et al 2004r) in th e southern hemisph ere, the 



extended Brightest Cluster Sample (eBCS) flEbehng et al 19981 : lEbeling et al 2000f ) in the 
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north, and the Clusters in the Zone of Avoidance (CIZA) (lEbehng et al 20021 : iKocevski et al 20071 ) 
sample along the Galactic plane. These are the most statistically complete X-ray selected 
cluster catalogs ever compiled in their respective regions of the sky. All three surveys are 
X-ray selected and X-ray flux limited using RASS data. The creation of the combined all-sky 



catalogue of 782 clusters is described in detail by (IKocevski et al 20061 ) and KA-BKE. 



We started with 3- year "foreground- cleaned" WMAP data (http://lambda.gsfc.nasa.gov) 
in each differencing assembly (DA) of the Q, V, and W bands. Each DA is analyzed sepa- 
rately giving us eight independent maps to process: Ql, Q2, VI, V2, Wl,..., W4. The CMB 
maps are pixelized with the HEALPix parameter A'"side=512 corresponding to pixels ~ 7' on 
the side or pixel area 4 x 10~^ sr. This resolution is much coarser than that of the X-ray 
data, which makes our analysis below insensitive to the specifics of the spatial distribution 
of the cluster gas, such as cooling flows, deviations from spherical symmetry, etc. In the 
filtered maps for each DA we select all WMAP pixels within the total area defined by the 
cluster X-ray emission, repeating this exercise for cluster subsamples populating cumulative 
redshift bins up to a fixed z. In order to eliminate the influence of Galactic emission and 
non-CMB radio sources, the CMB maps are subjected to standard WMAP masking. The 
results for the different masks are similar and agree well within their statistical uncertainties. 

The SZ effect (oc rig, the electron density) has larger extent than probed by X-rays (X- 
ray luminosity oc ng), which is confirmed by our TSZ study using the same cluster catalogue 
(AKKE). As shown in AKKE, KA-BKE contributions to the TSZ signal are detected out to 
5^30'. What is important in the prese nt context, is that t he X-ray emitting gas is distributed 



as expected from the ACDM profile (INavarro et al 19961 ) scaling as ne oc r ^ in outer parts. 



In orde r to be in hydrostatic eq uilibrium such gas must have temperature decreasing with 



radius (IKomatsu fc Seljak 20011 ). Indeed, the typical polytropic index for such gas would 
be 7 ~1.2, leading to the X-ray temperature decreasing as Tx oc rC~^ oc r~°-^ at outer 
radii. This Tx decrease a grees with simulations of cluster formation within the ACDM 



model (jBorgani et al 2004 ) and with the available data on the X-ray temperature profile 
(jPratt et al 20071 ). For such gas, the TSZ monopole (oc Txt) decreases faster than the KSZ 
component (oc r) when averaged over a progressively increasing cluster area. To account 
for this, we compute the dipole component of the final maps for a range of effective cluster 
sizes, namely [1,2,4, 6]6'x-ray and then the maximal cluster extent is set at 30' to avoid a few 
large clusters (eg. Coma) bias the dipole determination. We note that at the final extent our 
clusters effectively have the same angular radius of 0.5°. Our choice of the maximal extent 
is determined by the fact that the SZ signal is detectable in our sample out to that scale 
(AKKE), which is ~(3-4)Mpc at the mean redshift of the sample. Increasing the cluster 
radius further to 1° — 3°, causes the dipole to start decreasing with the increasing radius, as 
expected if the pixels outside the clusters are included diluting the KSZ signal (KA-BKE). 
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To estimate uncertainties in the signal only from the clusters, we use the rest of the map 
for the distribution and variance of the noise in the measured signal. We use two methods 
to preserve the geometry defined by the mask and the cluster distribution: 1) N^i central 
random pixels are selected outside the mask away from the cluster pixels adding pixels within 
each cluster's angular extent around these central random pixels, iteratively verifying that 
the selected areas do not fall within either the mask or any of the known clusters. 2) We 
also use a slightly modified version of the above procedure in order to test the effects of the 
anisotropy of the cluster catalog. There the cluster catalog is rotated randomly, ensuring 
that the overall geometry of the cluster catalog is accurately preserved. Both methods yield 
very similar uncertainties; for brevity, we present results obtained with the first method. 



2. Results 

Fig. [T] summarizes our results averaged over all eight DA's. We find a statistically 
significant dipole component produced by the cluster pixels for the spheres and shells ex- 
tending beyond z ~ 0.05. It persists as the monopole component vanishes and its statistical 
significance gets particularly high for the y-component. The signal appears only at the clus- 
ter positions and, hence, cannot originate from instrument noise, the CMB or the remaining 
Galactic foreground components, the contributions from which are given by the uncertainties 
evaluated from the rest of the CMB map pixels. The signal is restricted to the cluster pixels 
and thus must arise from the two components of the SZ effect, thermal and/or kinematic. 

The TSZ component, however, is given by the monopole term at the cluster positions 
and cannot be responsible for the detected signal. For the largest apertures it vanishes within 
the small, compared to the measured dipole, statistical uncertainty, and yet the dipole term 
remains large and statistically significant. This is the opposite of what one should expect 
if the dipole is produced by the TSZ component. Any random distribution, such as TSZ 
emissions, would generate dipole oc (tTx cos 9) which can never exceed (and must be much 
less than) the monopole component of that distribution, oc (tTx). On the other hand, any 
coherent bulk flow would produce dipole oc Vbuik(T cos^ 6*), which is bounded from below 
by the amplitude of the motion. Furthermore, we find significant dipole from (at least) 
-2mean=0.035 (135 clustcrs) all the way to 2;mean=0.11 (674 clusters); its parameters do not 
depend on the numbers of clusters, pixels used etc. Any dipole component arising from the 
TSZ term would depend on these parameters as it reflects the (random) dipole of the cluster 
sample and should thus decrease as more clusters are added in spheres out to progressively 
larger z. To verify this, we compute the expected monopole and dipole terms produced by the 
TSZ effect using the parameters of our cluster catalog as discussed in KA-BKE and recover 
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the monopole term fairly accurately when the /3-profile assumption is reasonable. The TSZ 
dipole component is then a small fraction of the monopole term. When normalized to the 
remaining monopole term in Fig. [T^ it is completely negligible compared to the measured 
dipole. Further, the TSZ dipole becomes progressively more negligible as more clusters are 
added in at higher 2;, and its direction varies randomly reflecting the random nature of the 
intrinsic cluster sample dipole on these large scales. All this is contrary to what we measure. 

We thus conclude that the dipole originates from the KSZ effect due to the bulk flow of 
the cluster sample. Our results indicate a statistically significant bulk-flow component in the 
final filtered maps for cluster samples in the z-bins from z <0.05 to <0.3 corresponding to 
median depth to z ~0.1, and it also persists when the dipole in shells is computed selecting 
only clusters at z >0.12 (the median redshift for this sub-sample is ~0.18). Fig.[T^ shows that 
the bulk flow results in a CMB dipole with little variation - within the statistical uncertainties 
- between ^median — 0.03 and iJ0.12. The monopole component reflects the residual TSZ 
contribution which is very small for the maximal cluster aperture as Fig. [1^ shows. (At 
lower z there may still be some residual TSZ component, which would be consistent with 
the more nearby clusters having a larger angular SZ extent than the more distant ones). 

To translate the CMB dipole in /iK into Vbuik in km/sec, we generated CMB tempera- 
tures produced by the KSZ effect for each cluster and estimate the dipole amplitude, Ci.ioo, 
contributed by each 100 km/sec of bulk-flow (KA-BKE). The results are shown in the last 
column of Table 1 of KA-BKE for the central values of the direction of the measured flow; 
within the uncertainties of (/, 6) they change by at most a few percent. A bulk flow of 100 
km/sec leads to \/C\^\qq ~ l/iK for unfiltered clusters assuming the /3-model; this corre- 
sponds to an average optical depth of our cluster sample of (r) ~ 10~^ expected for a typical 
galaxy cluster. For NFW clusters the value of Ci.ioo would be smaller. Filtering reduces the 
effective r by a factor of 3. Since a /3-model provides a poor fit to the measured TSZ com- 
ponent outside the estimated values of 6'x-ray, we compute Ci^ioo within that aperture where 
the central value of the bulk-flow dipole has approximately the same value as at the final 
apertures. Due to the large size of our cluster sample {Nd ~130-675), the random uncertain- 
ties in the estimated values of Ci^ioo should be small, but we cannot exclude a systematic 
offset related to selection biases affecting our cluster catalog at high z. Such offset, if present, 
will become quantifiable with the next version of our X-ray cluster catalog (in preparation) 
using the empirically established SZ profile rather than the current (3-mode\ to parameterize 
the cluster TSZ. The good agreement between the various TSZ-related quantities shown in 
KA-BKE for 9sz=0x-Tay and the observed values for both unfiltered and filtered maps sug- 
gests, however, that these systematic uncertainties are likely to be small. They only affect 
the accuracy of the determination of the amplitude of the bulk flow, but not its existence 
established by the CMB dipole at the cluster locations. Since the filtering removes r in the 
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outskirts of clusters more effectively, a larger amount of power is removed in the /3-model 
when the cluster SZ extent is increased beyond 6'x-ray5 than in the steeper profile measured 
by AKKE. Thus the effective r is possibly underestimated by using a /3-model, but it cannot 
exceed (and must be much less than) the calibration obtained from the unfiltered /5-model. 

3. Cosmological implications 

Conventionally, the entire peculiar velocity field is assumed to be driven by the peculiar 
gravitational potential. For a given cosmological model, the details of the velocity field also 
depend on the window function of the dataset. Constructing the precise window function is 
beyond the scope of this paper, but the overall conclusions would be insensitive to its details 
because the amplitude and coherence length of the measured flows are quite unexpected 
within the concordance ACDM model. Fig. [T]F shows the rms prediction, ay, of the concor- 
dance ACDM model. If produced by gravitational instability within the concordance ACDM 
model, the motion would require the local Universe out to ~ 300/i^^Mpc to be atypical at 
the level of many standard deviations of the model. Indeed a value of a/C%ioo ~ 3/iK is 
required to reach peculiar velocities of order 100 km/sec on the relevant scales. This is much 
greater than Ci^ioo deduced from the unfiltered y^-i ay data and even then it would be difficult 
to explain the approximate constancy of the measured dipole with depth. 

Cosmic variance does not change these conclusions significantly. For a Gaussian density 
field the peculiar velocity distribution on linear scales is Maxwellian, with the probability 
density of measuring a 1-D bulk velocity p{y)dV o^V'^ e^^{-l.bV^ / al)dV . The probability 
of finding a region with V <Vq is then P(Vo)=r(|, where F is the incomplete gamma- 
function normalized to F(?t,, cxd)=1. The (68%, 95%) c.l. require Vo=(1.08,1.6)(Ty; the shaded 
area in Fig. [T]F shows the 95% c.l. region. In the concordance ACDM model (Jv/=(150,109) 
km/sec at (200, 300)/i~^Mpc, so 95% of cosmic observers should measure bulk flow velocities 
less than (240,180) km/sec at these scales. To make these numbers consistent with our 
measurements - at these scales alone - would require ■\/Ci,ioo~2/iK. This is much higher 
than even the calibration values for unfiltered data (computed for the /5-model which leads 
to a larger Ci^ioo value than the NFW-profile clusters) and cannot be accounted for by any 
systematic uncertainties in our calibration procedure. Indeed, the distribution of errors is 
approximately Gaussian (see Fig. 5 of KA-BKE), so the probability of measuring velocity 
Vi ± ej at scale is Vi oc PiV) e'x.\)[—^^^^^^\dV (e.g. Gorski 1991). The probability 
of several such independent measurements is the product of Pj's; for the numbers plotted 
in Fig. If this probability is completely negligible. (The measurements in Fig. If are not 
strictly independent since each subsequent z-hm. also contains the clusters from the previous 
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bin; nevertheless each Vi is so small that the overall probability is still negligible). 

The coherence length of the measured bulk flow shows no signs of convergence out to 
^300/i~^Mpc, and it is quite possible that it extends to much larger scales, possibly all 
the way across our horizon. An interesting, if exotic, explanation for such a "dark flow" 
would come naturally within certain inflationary models. In general, within these models 
the observable Universe represents part of a homogeneous inflated region embedded in an 
inhomogeneous space-time. On scales much larger than the Hubble radius, pre-inflationary 
remnants can induce tilt inc luding CMS anisotropies generated by the Grischuk-Zeldovich 



(iGrischuk fc Zeldovich 19781 ) effect (ITurner 199ll : iKashlinsky et al 19941 ). These can arise 



from the parts of space-time that inflated at different times and rates and would manifest 
themselves mainly in the quadrupole component, Q: an inhomogeneity of amplitude ^x, ~ 1 
at a distance L ^ cHq^ generates a quadrupole Q ~ Sl{cHq^ / . Consistency with the 
observed low value of Q would require a sufficiently large number of the inflation's e-foldings, 
making the Universe flat to within |1 — fitotail < Q and causing t he scale of inhomogene - 



ity L to become comparable to the curvature radius (> 500cHq^) (IKashlinsky et al 19941 ). 
Such a tilted universe would lead to a uniform flow across the observed horizon due to the 
density gradient produced by this superhorizon mode. The bulk motion would have an 
ampli tude of f ~ cSr.jcHn^ / L) and would not generate a primordial dipole CMB compo- 



nent (iTurner 19911 ). Since the quadrupole produced by such a pre-inflationary remnant is 
Q ~ {v/c){cHq^/L), it is possible for such inhomogeneities to generate the required motions 
and be consistent with the observed value of Q and flatness. Although it would require acci- 
dental alignment, the contribution from such an inhomogeneity to CMB anisotropies via the 
GZ effect might, interestingly, also explain the observed low value of the CMB quadrupole 
(and possibly also octupole) compared to the concordance ACDM model. This explanation 
for the measured bulk flow would, however, still require peculiar velocities generated by grav- 
itational instability acting on the ACDM density fleld, which would provide a random compo- 
nent around the uniform bulk flow. On sufficiently large scales, such a flow would be "cold " 



in the sense that it would be characterized by a large Mach number (jOstriker &: Suto 19901) 



which may be measurable in future cluster surveys (lAtrio-Barandela et al 20041 ): the Mach 
number then should increase linearly with scale on scales ^ 100/;.~^Mpc. On smaller scales, 
there may be non-negligible contributions to the flow from peculiar motions generated by 
the gravitational instability caused by local matter inhomogeneities. This can lead to a 
non-alignment with the flow at lower z in general agreement with the trends in Fig. 1. 
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Fig. 1. — Upper panels: (a)-(d) Monopole and dipole terms for 9sz = min[(l, 2, 4, 6) x 
^x-ray, 30'] with 1-(T Standard deviations. Values at the maximal aperture, which have the 
lowest monopole term, are marked with filled circles. The averaged monopole and dipole 
components are weighted with statistical uncertainties. The statistical significance of the 
KSZ component improves as more of the cluster pixels producing the signal arc included at 
higher z. The noise of our measurement of the dipole at 1.8(Ai"ci/100)~^/^/xK with three- year 
WMAP data is in good agreement with the expectations of KA-B. Lower panel, (e) - Outer 
2;-bins with signal measured at ^2a. Filled circles show the values from Table 1 of KA-BKE 
at the maximal aperture vs the median z; open symbols show the same vs the mean z. The 
two symbols are connected to show the uncertainty in the scale on which the fiow is probed. 
Signal recovered at ^sz = min[(l,2,4) x ^x-ray, 30'] is shown with la error bars; from left 
to right in order of increasing aperture. The values are slightly displaced around the true 
-^median for clearer display, (f) - Comparison between theoretically expected bulk fiow and 
the measurements. The rms bulk velocity for the concordance ACDM model which best 
fits the WMAP 3-year data for top-hat (solid line) and Gaussian (dashes) windows; shaded 
region marks the 95% cl from cosmic variance. The results of this study, translated into 
km/sec using -y/ Ci,ioo = 0.3/iK, arc shown with 1-a errors vs the mean/median rcdshift of 
the clusters in each cumulative z-hin. The horizontal bars connect Zmenn with 2;median- The 
results in shells from Table 1 in KA-BKE are omitted in this comparison because of the 
theoretical windows plotted, but they show that the motion extends to mean redshift ^0.18 
well beyond the horizontal range of the figure. 



